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ABSTRACT 

We report on the results of HST optical and UV imaging, Spitzer mid-IR photometry and optical spectroscopy 
of a sample of 30 low-redshift (z^O.l to 0.3) galaxies chosen from the SDSS and GALEX surveys to be 
accurate local analogs of the high-redshift Lyman Break Galaxies (LBGs). The Lyman Break Analogs (LBAs) 
are similar in stellar mass, metallicity, dust extinction, SFR, physical size and gas velocity dispersion, thus 
enabling a detailed investigation of many processes that are important in star forming galaxies at high redshift. 
The main optical emission line properties of LBAs, including evidence for outflows, are also similar to those 
typically found at high redshift. This indicates that the conditions in their interstellar medium are comparable. 
In the UV, LBAs are characterized by complexes of massive clumps of star-formation, while in the optical 
they most often show evidence for (post-)mergers and interactions. In six cases, we find a single extremely 
massive (up to several x 10^ A/0) compact (radius ~ 10^ pc) dominant central object (DCO). The DCOs 
are preferentially found in LBAs with the highest mid-IR luminosities (^24/^™ = 10^" '^ to 10^^'^ Lq) and 
correspondingly high star formation rates (15 to 100 Mq per year). We show that the massive star-forming 
clumps (including the DCOs) have masses much larger than the nuclear super star clusters seen in normal late 
type galaxies. However, the DCOs do have masses, sizes, and densities similar to the excess-light/central- 
cusps seen in typical elliptical galaxies with masses similar to the LBA galaxies. We suggest that the DCOs 
form in the present-day examples of the dissipative mergers at high-redshift that are believed to have produced 
the central-cusps in local ellipticals (consistent with the disturbed optical morphologies of the LBAs). More 
generally, the properties of the LBAs are consistent with the idea that instabilities in a gas-rich disk lead to very 
massive star-forming clumps that eventually coalesce to form a spheroid. Finally we comment on the apparent 
lack of energetically significant AGN in the DCOs. We speculate that the DCOs are too young at present to 
be growing a supermassive black hole because they are still in a supernova-dominated outflow phase (age less 
than 50 Myr). 

Subject headings: galaxies: starburst — galaxies: active — galaxies: bulges — galaxies: peculiar — galaxies: 
high redshift 



1. INTRODUCTION 

One of the most remarkable features of high redshift (z ^ 
2 — 6) galaxies is their great potential for vigorous star for- 
mation. Moreover, observations with the Hubble Space Tele- 
scope (HST) indicate that that star formation largely occurs in 
an ext remely "clumpy" and compact mode (e.g. Cowie ( 
T99I iFranx et alJ \\Wk rElmegreen et alJ 12005; .Stark 
200a). These clumps are reminiscent of star-forming com 
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plexes in local HII regions, but their total masses and scale 
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sizes are several orders of magnitudes larger As shown by 
lElmegreen et al. (2009), clumps with masses of 10^^^ Af© 
and sizes of >1 kpc are a very typical constituent of ex- 
tended , irregular galaxies at z — 1 — 4 in the Ultra Deep Field 
(UDF, i Beckwith et al.ll2006l) . The relatively high gas veloc- 
ity dispersions observed in some high redshift galaxies (e.g. 
Eaw et al. 2007, 2009; Forster Schreiber et al. 2009) coupled 
with their large inferred gas column densities and gas frac- 
tions can explain the formation of the clumps and their bulk 
properties from gravitational instabilities in gas-rich, Jeans- 
unstable disks. This process may be instigated, for example, 
by the continuous accretion of gas o r by the dehvery of gas- 
rich material in small merger events ( iBournau d & ElmegreenI 
2009; Elmegreen et al. 2009). In any case, the compactness of 
the stellar clumps furthermore implies that they can form rel- 
atively unhampered by velocity shear or supernova feedback. 
Simulations show that for some galaxies the coalescence of 
clump systems may lie at the root of the formation of galaxy 
bulges, and possibly even their supermassive nuclear black 
holes if it is assumed that each clump gave rise to its own 
interme diate mass black hole prior to their migration to the 
center (|Noguchilll99l Ilmmelietal.ll200l; lElmegreen et"an 
l2008allb ir^ 

It would be extremely valuable to study similar objects in 
nearby galaxies where multi-waveband data with relatively 
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high spatial resolution and sensitivity would allow the rele- 
vant physical processes to be probed in much greater detail. 
While local starburst galaxies used in previous studies have 
some similarities to star-forming galaxies at high redshift, 
there are some important differences: local starbursts hav- 
ing SFRs of >100 AIq yr^^ are very dusty systems (the so- 
called Ultra-Luminous Infrared Galaxies, ULIRGs). The in- 
tense star formation occurs in compact, sub-kpc scale regions 
in the centers of interacting g alaxies, and are usually incon- 
spicuous m the far-UV ( Heckman et al.|[T99l iGoTdader et alj 
12002). Although very dusty starbursts clearly make u p a non 



negligible fraction of high redshift galaxies (e .g. Hughe s et al 
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iHuang et al.l l2005t IPapovich et al 
purgarella et al. 2009; Daddi et al 



20091). most of the 



systems contributing to the cosmic SFR at high redshift are 
accounted for in deep UV-selected sa mple s of Lyman Break 
Galaxies (LBGs) (e.g. Redd v et all I2006t iReddv & Steidell 
|2009; Bouwens et al. J009), indicating that the amount of 
dust extinction in LBGs is typically much lower: the average 
ratio of their bolometric IR-to-UV luminosities, Lm/ Ljjv, is 
much lower (~ 10 x ) compared to local galaxies of the same 
Lboi (e.g. B uat et al.l2007l: IReddv et alJ2 006^iBur gareUa et alj 
12007). Locally calibrated dust corrections that are based 
on the reddenin g of th e UV continuum slope as derived 
iMeure r et al.l 
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19991) seem to apply to ty pical LBGs 



^eddv et al. 2006 ) with some exceptions (e.g. iBaker et alj 
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to show that, despite their UV selection, LBGs at z ~ 2 — 4 
still put out most of their energy in the mid- and far-IR 
dBouwens et al.ll2009l) . 

The star formation in LBGs typically occurs over the entire 
extent of the galaxy ranging from a few hundred p c to several 
kpc (e.g. [F ranx et al. 1997; Giavalisco 2002; Bo uwens et alj 
l2004 lFerguson et al. 2004 iStark et al.i ,2008). The relatively 
low extinctions and high UV surface brightnesses seen in 
LBGs are characteristic of local blue compact dwarf galax- 
ies (BCDs). However, typical BCDs have star formation rates 
up to two orders of magnitude small er than LBGs and occur in 
galaxies wi th much lower mass (e.g. lTelles & Terlevichll997t 
iHopkins eT al. 2002), although some may resemble fainter 
LBGs at high redshift in their basic properties (see, e.g., 
lOstUn et al.i.200 9). The main differences between LBGs and 
typical local starbursts studied previously may be due to a sys- 
tematically lower metaUicity and higher gas-mass fraction in 
the LBGs and/or d ue to a different t riggering mechanism (e.g. 
lErb et alJ l2006a"F: 'L aw et al.ll2009h . 

The "Lyman break analogs" (LBA) project was designed 
in order to search for local starburst galaxies that share typ- 
ical characteristics of high redshift LBGs jHeckman et alj 
2005). The LB A sample and main prope r ties are given in 
Heckman et al. (2005), Hoopes et al. ( 2007b. iBasu-Zvch et alj 
P007, 2009) and Overzier et al. (2008, Paper I), and we note 
that the LBAs are identical to the sample of "supercompact 
UV-luminous galaxies" (ScUVLGs) referred to in the papers 
listed above. In brief, the UV imaging survey performed by 
the Galaxy Evolution Explorer (GALEX) was used in order 
to select the most luminous (Lpuv > 10^°'^ Lq) and most 
compact {Ipuv > 10^ Lq kpc^^) star-forming galaxies at 
z < 0.3. Galaxies having these high Lpuv are already ex- 
tremely rare (local space density of ^ 10^^ Mpc~^), and 
the surface brightness requirement further reduces their ab- 
solute number density by a factor of ~ 6. These objects 
tend to be much more luminous in the UV than typical local 



starburst galaxies studied previously dHeckman et alj [T998I: 
iMeurer et alll 19971 119991) consistent with high SFRs and rel- 
atively little dust extinction. The median absolute UV mag- 
nitude of the sample is -20.3, corresponding to ~ Q.5L*^^, 
where L't^^ is the characteristic luminosity of LBGs at z ^ 3 
(Steideletal. 1999, i.e., Muqo^ab = -21.07). Analysis of 
their spectra from the Sloan Digital Sky Survey (SDSS) and 
their spectral energy distributions from SDSS and GALEX 
subsequently showed that the LBAs are similar to LBGs in 
their basic global properties, including: stellar mass, metal- 
licity, dust ex tinctio n, SFR, physical size and gas velocity dis- 
persion. In Overzier et aXI (12008 ) we showed that most of the 
UV emission in a preliminary sample of 8 LBAs originates in 
highly compact burst regions in small, clumpy galaxies that 
are also morphologically similar to LBGs. We demonstrated 
that if LBGs at high redshift are also small merging galaxies 
similar to the LBAs, unfortunately this would be very hard to 
detect given the much poorer physical resolution and sensitiv- 
ity. Furthermore, for the present paper it is important to point 
out that the LBAs also occupy an "offset" region in the main 
optical emission line diagnostics dia gram (i.e. log([NII1^ a) 
versus log([OIII]/H/3), or the "BPT" ffialdwin et al.lll98 i') di- 
agram) analogous to galaxies at high redshift. This suggests 
that similar physical conditions may apply to their interstellar 
medium (ISM) as well. In summary, the LBAs appear indeed 
good local analogs to the LBGs. 

In this paper we will present new results from our HST and 
Spitzer imaging campaigns and spectroscopic follow-up in or- 
der to highlight some of the remarkable properties of the LBA 
sample: we will compare the UV-optical sizes and morpholo- 
gies, showing evidence for massive star-forming "clumps" 
that dominate the UV (§3). We explore possible connections 
between the starburst regions and the general properties of the 
optical emission line gas (§3.5 and §4), as well the relation- 
ship between the formation and evolution of massive stellar 
clumps in the context of the main structural components of 
young galaxies (§5). We conclude the paper with a summary 
and final remarks (§6). Forthcoming work on the LBA sam- 
ple currently in preparation includes a detailed morphologi- 
cal analysis, multi-waveband modeling, mid-IR spectroscopy, 
VLBI radio imaging. X-ray spectroscopy, and integral field 
and long-sht optical and near-IR spectroscopy. 

2. DATA 
2.1. HST and Spitzer Imaging 

We observed 30 LBAs using HST (Programs 10920/1 1 107, 
PI: Heckman). In order to probe morphologies in the UV, 
twenty-four objects were observed using the ACS Solar Blind 
Channel (SBC) in the filter F150LP (A^ « 1614A), six ob- 
jects were observed using the High Resolution Camera (HRC) 
in the filter F330W (Ac ~ 3334A). The exposure times per 
object were typically ^2500 s. These UV images offer a 
tremendous gain in resolution over that of the GALEX im- 
ages («G'.'07 compared to «5" for GALEX). Matching obser- 
vations in the rest-frame optical were carried out either using 
the Wide Field and Planetary Camera 2 for ^^3600 s in F606W 
(Ac ~ 6OOIA; twenty-four objects) or using the ACS Wide 
Field Channel for -2200 s in F850LP (Ac « 917GA; seven 
objects). The ACS and WFPC2 observations were divided 
into, respectively, three and six dit hered exposures and we re 
combined using multidrizzle dKoekemoer et al.ll2002h . 

We make use of photometric data obtained with the Multi- 
band Imaging Photometer (MIPS) on the Spitzer Space Tele- 
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scope at 24 /im (Program ID 20390, PI: C. Hoopes). Photom- 
etry was obtained by performing point source extractions on 
the post-Basic Calibration Data. 

2.2. SDSS and VLT/FLAMES Spectroscopy 

We further use integrated optical emission line fluxes and 
flux ratios measured from the SDSS spectra. In order to ac- 
count for stellar absorption lines, the emission-line parame- 
ters were obtained from Gaussian fits after subtracting the 
continuum based on the best-jitting spectrum using stellar 
population models (see lTremonti et al] i2004). 

The fraction of the total (line and continuum) light that is 
included within the 3" diameter SDSS fiber apertures is found 
to be 0.6 ± 0.1, estimated from the difference in fluxes mea- 
sured in a 3" aperture compared to the total SDSS photomet- 
ric aperture'^. We estimate that this estimate is in fact a lower 
limit to the fraction of UV and line emission included within 
the SDSS fiber aperture, based on the compactness of the star- 
forming regions in the HST images and follow-up spectra. 

In Sect. 14.11 we present a preview of data" taken 
with the Fibre Large Array Multi Element Spectrograph 
(FLAMES) on the European Southern Observatory (ESO) 
Very Large Telescope (VLT) during the night of Oct 30 2008. 
We used the GIRAFFE/ARGUS integral field spectrograph 
(IFS) to obtain seeing-limited spectra resolving the Ha and 
[NII]AA6548A,6584A line complex at a spectral resolution of 
R ~ 9000 and a spatial plate scale of 0'.'3 or 0'.'52 pixer^ 
(depending on the seeing). Each source was observed in 4 
exposures of 400 s each. The data was pipeline reduced and 
further cleaned of cosmic rays and the background using our 
own routines. 

More details on the LBA sample selection, observations 
and all measur ements used in this paper are given else- 
where (Paper I: iHoopes et alj|2007t iBasu-Zvch et"aLl 120071: 
lOverzier et al.ll20'09h . The object coordinates, redshifts and 
basic quantities used in this paper are listed in Table[Tl Emis- 
sion line measurements are listed in Table |2l 

2.3. Measuring Stellar Masses 

We will use stellar mass estimates that were obtained from 
the SDSS/DR7 stellar mass catalogue'^. These masses are 
derived from fits to the SDSS photometry u sing a lar g e grid 
of models constructed using the Bruzual & CharlotI (|2003l 
BC03) stellar synthesis library. Although the BC03 models 
may not include the most up-to-date treatment of the ther- 
mally pulsating (TP)-AGB phase, it has been shown that this 
phase mostly affects the stellar mass-to-light ratio in the rest- 
frame near-infared H and K bands (e.g., see Eminian et al] 
l2008i) . Because near-infared data is not being used in the fits, 
we believe that our stellar mass estimates are relatively accu- 
rate. However, in order to further reduce the uncertainty in 
these galaxy stellar mass estimates, we will derive a "best" 
mass by taking the average of the SDSS photometric mass 
and a dynamical mass that we calculate from the Ha emis- 
sion line velocity dispersion and the optical half-light radius. 
The two mass estimates typically agree to within ^0.3 dex 
with no systematic differences (see Table[T]i. 

The fraction of fiber-to-total-light in the SDSS r'-band is estimated by 

f fiber = 10-'^-*'^"'r;fib.r--m,..P,tro) 

" ESO Program 082.B-05 12(A) 

Available at http : / /ww w . mpa-garching . mpg . de/SDSS/DR7/| 



2.4. Measuring Star Formation Rates 

We will make use of three different estimates of the inte- 
grated star formation rate. We use the most recent c alibrations 
from literature as tab ulated in lCalzetti et al.l (12009b . and valid 
for a lKroupal ( l2008l) initial mass function (IMF). The result- 
ing SFRs are lower by a factor of ~ 1 .5 compared to a Salpeter 
(1955) IMF 

2.4.1. Ha 

SFR/fcf is calculated from the Ha line luminosity and ap- 
plying a correction for dust based on the Balmer decrement 
XfHa/fHfj) following the recipe in (.Calzettii,2001l) : 

SFRff„,o [Mq yr-i] = 5.3 x IQ-^^^h^.o [erg s-^], (1) 

where Lh^q = ^Ha,ofcslO° '''^"°, and 

AHo.^KHa)E{B~V)gas, (2) 
_ logio[(fea//g/3)/2.87] 
^ - o.4[fc(i/a) - fc(iJ/3)] ' 

with k{Ha) = 2.468 and k{Ha) - k{H(j) = 1.17. We 
apply a small correction factor of typically ~ 1.7 due to the 
flux expected outside the SDSS fiber. 

2.4.2. Ha+24 fim 

An improvement to the above method is given by 
SFR^fa-|_24, which is based on a combination of the uncor- 
rected Ha luminosity (Lffa,obs), ™d the 24/im luminosity 
(L24 = ^^24^24) related to the emission of dust heated by 
young stars. It has been shown that SFRHa+24 is a robust in- 
dicator for the total SFR calibrated a gainst Paa measurem ents 
in a large sample of nearby galaxies dCalzetti et al.ll2007h : 

SFRh„+24 [Mq yr-1] = (3) 
5.3 X 10-4' X (Lh„,o6s +O.O3IL24) [crgs-i]. 

As the exact, individual infrared spectral energy distribu- 
tions are currently uncertain, we choose not to apply any K- 
corrections to the observed 24/im fluxes. Based on the av- 
erage corrections found for a collection of spectral templates 
presented by Rieke et al. (2009, see their Figure 9) we esti- 
mate that the actual SFRs could be between 0.0 and 0.3 dex 
higher depending on the redshift, spectral shape, and total IR 
luminosity. 

24.3. FUV 

We will make use of SFRi?t/v.o> which is based on the 
dust-corrected far-UV luminosity using the empirical cor- 
relation between the attenuation of the UV and the IR 
emission from dust heated by obscured stars (the "(3-IRX 
relation", iMeurer et al.l Il995h . The correlation has been 
shown to work well for relatively unobscured star-forming 
and starburst galaxies at both low and high r edshift (e.g. 
Meure retal. 19 95, 1999; ISeibert et al.1 l2005t iReddv et alj 
I2006t ISalim et al. 2007), but it must be noted that the situa- 
tion is more complex for galaxies that have more complicated 
star f ormation histories (e.g. Kong et al. 2004; Johnson et all 
(2007!), very young ages (e.g. iReddv et al.. ,200fe .Siana et alj 
2008 ) or that are heavily obscured (e.g T lMeurer et al. 
IGoldader et al.ll2002l: iReddv et al.l l2006'). We use the calibra- 
tion between UV colour and Apuv from .Trever et al.. (,2007i) . 
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Fig. 1. — False-color HST images of the LBA sample showing the (rest-frame) UV in blue/purple and optical in yellow/red. The images measure 6"x6". 
The UV images were rebinned (2x2) to match the pixel scale of the optical images, and convolved with a Gaussian kernal with ff/l (FWHM). Although most 
objects are highly compact in both the UV and optical, a small subset consists of a very bright unresolved component in the middle of an extended, low surface 
brightness disk. The images demonstrate a wide range of complex morphologies often suggestive of interactions and (post-)merging. See text for details. 



Afuv = 4:.05{inFuv 
SFRfuv [Mq yr-1] 



rriNuv) — 0.18 mag, finding: 
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Using the spectral slope f3 = 2.321ogio("^i^(7V — rriNuv) — 
2.0 with f\ oc A'' we further applied a small ii'-correction 
to obtain the monochromatic luminosity at rest-frame 1500A, 

1500A 

3. RESULTS 

3.1. UV-optical morphologies 

In Fig. [T]we show 6"x6" UV-optical false-color postage 
stamps of 30 LBAs, including the 8 objects studied in Pa- 
per I. In all cases, both the UV and the optical light is dom- 
inated by the intense light from young massive stars in very 
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Fig. 2. — Stellar mass versus the half-light radius measured in the UV (left 
panel) and optical {right panel) images. Sohd sq uares mark the subset of 
DCOs, in anticipation of results obtained in Section U!2l 

compact starburst regions. The optical images furthermore 
allow us to see the faint emission from older stars display- 



Local Analogs to Lyman Break Galaxies 



5 



032845 
1 


04ffi?oe 


082001 

♦ 


082550 

/ 


083803 
1 


^^^^ 




Fig. 3. — Illustration of HST optical morphologies of typical LBAs having 
several star-forming clumps of comparable brightness. These are to be con- 
trasted with the six LBAs having a dominant compact central object shown 
in Fig.|4] The images measure 10" X 10". 
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Fig. 4. — HST optical morphologies of the six LBAs having a dominant 
compact central object. The images measure 10" X 10". 

ing a wide range of (often complex) morphologies, which in- 
clude clumpy or chain-like galaxies, interacting pairs, a ring 
galaxy, and perturbed (spiral) disks. In Fig. |2]we plot the 
half-light radii versus the stellar masses. The global median 
half-light radii of the galaxies are «1 kpc and «L5 kpc in 
the UV and optical respectively. The galaxy mass estimates 
range from a few times 10^ to < 10^^ M© (Table [T]i. The 
sizes and the stellar masses of these galaxies correspo nd well 
to those of typical Lyman Break Galaxies (e. g. Shapl ev et al. 
2001'; 'Papovich et al."2001t lGiavaliscoll200l iFerguson et al. 
2004; Stark et al. 2009). 

As noted already in Paper I, a significant fraction of both 
the UV and optical emission is seen to come from luminous, 
very compact clumps. The much larger size of the sample 
presented here allows us to investigate this in more detail. In 
Figs. [3] and |4] we show some close-ups of the rest-frame op- 
tical HST images at high contrast. LBAs typically consist 
of a collection of centrally-located compact clumps having a 
range in magnitudes (e.g.. Fig. [3]l. However, in a minority 
of cases the LBA contains a single, dominant luminous point- 
like source at or near the center of the galaxy (i.e., 005439, 
021348, 080232, 080844, 092159, and 210358). In these 
sources, the central source typically contains ~50% (20 and 



77% being the lowest and the highest, respectively) of all the 
UV light. In the optical, the contribution from the host galax- 
ies is larger, somewhat reducing the light fraction of the cen- 
tral source to 20 — 30%. Here we have used an aperture 
of 0'.'2 radius for the central point-like source and an aperture 
of 10" to estimate the total light of the galaxy. Because of 
their visual prominence and in anticipation of some of their re- 
markable properties that will be investigated in the following 
sections, we will henceforth refer to this subset of the clumps 
as "dominant central objects" (DCOs). A collage of images 
of the six objects hosting DCOs is shown in Fig. |4] 

3.2. Color-magnitude diagrams: estimated clump masses 

We have measured the magnitudes and colors of each of 
the clumps and DCOs identified by eye inside circular 0'.'4 di- 
ameter apertures. The backgrounds were determined inside 
a local annulus and subtracted. We did not make a correc- 
tion for the (small) amount of flux expected to fall outside 
the aperture due to the HST PSF In Fig. |5] we show the 
color-magnitude diagrams for all objects for which two fil- 
ters were available, either FUV and V (shown on the left) or 
U and z (shown on the right). Open and filled squares indicate 
the measured values, while the dotted lines indicate the dust- 
corrected values based on the Balmer decrement {fHalfHp) 
using a constant relation between reddening of the gas and the 
stellar continuum'^ ( ICalzettil200Ih . Unfortunately we need to 
assume that the Balmer decrement as measured from the in- 
tegrated SDSS spectrum applies to the much smaller clumps 
as well, although this is not necessarily correct. The results 
are tabulated in Table [3] The clumps are very luminous with 
a median absolute magnitude in V (approximately rest-frame 
B) of ~ -20 mag, and very blue {FUV -V ~ 0.0). The 
DCOs are the most luminous having Mb ^ -20 to -22 mag. 
The small gap in the magnitude histograms in the top panels 
of Fig. |5]may suggest that the DCOs are not just a bright tail 
of the general clumps distribution, but perhaps a distinct class 
of object. 

A c omparison with Starburst99 5.1 (IVazquez & Leithered 
I2005h stellar synthesis models indicates that the stellar masses 
of the clumps shown in Fig. |5]are typically a few times 10^ 
to a few times 10^ Af©, while the DCOs have typical masses 
near or in excess of 10^ Af© (see Table |3]l. It is important to 
note that these mass estimates are relatively insensitive to dust 
extinction, since the dust extinction vector runs nearly paral- 
lel to the age vector in these plots. As shown in Fig. |6] the 
masses are relatively insensitive to the star formation history 
(SFH) that we adopt. The figure also illustrates that it is much 
less straightforward to determine the typical stellar age of 
the clumps without accurate information about the extinction 
and SFH. For our adopted extinction corrections, the instan- 
taneous burst models imply ages in the range from ^5 to 100 
Myr (red points in Fig. |6]l. Alternatively, if we adopt a contin- 
uous SFH the ages would be substantially higher, in the range 
--100-1000 Myr (blue points), while the SFRs would have 
to be on the order of '--^10 M© yr^^ in order to explain the 
high luminosity of the brightest clumps (left blue track in Fig. 
|5j. The masses would be ~ 2x higher on average compared 
to those derived for instantaneous bursts (Table [3]l. However, 
there are several arguments to make in favor of an instanta- 
neous SFH (or a SFH having a very short e-folding time). We 
note that at least three of the sources show prominent Wolf- 
Rayet features indicating the presence of very massive stars 

'3 E(B - V), = 0.44E(B - V)gas 



6 



Overzier, R.A. et al. 




Fig. 5. — Bottom panel: Color-magnitude diagram of star-forming clumps identified in the HST WFPC2 images. Solid points show measured values, while 
dotted lines connecting to open points indic ate dust-co rrected values based on the Balmer decrement and assuming a constant relation between reddening of the 
gas and the stellar continuum following C alzettil 420011) . Coloured lines indicate population synthesis tracks modeled using Starburst99 for either constant bursts 
with a SFR of 1 and 10 A/q yr~^ (right and left blue tracks) or an instantaneous burst normalized to a total stellar mass of 10* Mq and 10^ Mq (right and left 
red tracks). Ages in Myr have been indicated along the tracks. Blue dashed lines indicate the ages at which the constant burst models reach a total stellar mass 
of 10* Mq and 10^ Mq. The models used have a LMC metallicity, include the nebular continuum, and were redshifted to z = 0.15. The red arrow indicates 
the reddening vector corresponding to an extinction of Ay = 1.0 mag. The six DCOs are indicated by the filled squares. Top panel: Histogram of the absolute 
magnitude distribution corresponding to the clumps shown in the bottom panel. The measured (dust-corrected) distribution is indicated by the shaded (open) 
histogram. Typical massive star clusters in nearby galaxies have My — —(16 — 14) that would be as bright as My ~ —(19 — 17) when observed at the young 
age of ~10 Myr. The LEA clumps, and the DCOs in particular, are generally much brighter than typical massive star clusters. See Sect. [5]for discussion. 
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Fig. 6. — Clump masses and ages estimated from the color-magnitude di- 
agrams presented in Fig. [3] assuming either an instantaneous star formation 
history (red points) or a continuous star formation history (blue points). The 
masses estimated using the continuous model are only 2 X higher on av- 
erage compared to those derived for instantaneous bursts, while the ages es- 
timated using the continuous model are significantly higher compared to the 
instantaneous model. 



at a special phase only Myr after the onset of a massive 



burst ("e.g. lBrinchmann et al.ll2008bh . This would favor the in- 
stantaneous burst model. Also, we note that the very compact 
sizes of these objects (and thus the short dynamical times) are 
consistent with a short formation time (tform ^ tdyn)- 

3.3. Clump sizes 

It is clear from simple visual inspection of the HST images 
that the DCOs are marginally-resolved at best (note the strong 
diffraction spikes associated with the DCOs in Fig. 3). Many 
of the other clumps appear either as faint point-like sources or 
larger spatially resolved regions (see examples in Fig. |3]l. 

To obtain the best constraints on the sizes of the DCOs, we 
have first examined the radial flux profile of the three bright- 
est examples, using the ACS UV images as they have smaller 
pixels and a sharper PSF compared to the optical images. The 
result is shown in Fig. [T] Each panel shows the measured 
count rate (plusses) compared to the measured profile of a 
mode l PSF (dash ed line) modeled using the Tiny Tim 3.0 soft- 
ware (' Kristi[l995! ). We also simulated a range of 2D Gaussian 
models having a FWHM of 0'.' 005-0'.' 125 and convolved them 
with the model PSF. The results are shown using coloured 
lines. In all three cases, the data is most consistent with an 
object no larger than « 0'.'075 (FWHM), corresponding to ef- 
fective radii no larger than ~70-160 pc at z = 0.1 — 0.3). 

In order to be consistent across our entire clump sample 
and with our earlier definition of the clumps, we have mea- 
sured the physical sizes of the brightest clumps in each LBA 
by measuring the optical half-light radii of each clump within 
the same fixed angular apertures used in Fig. |5] The radii were 
deconvolved assuming a PSF of O'.'l (FWHM). In Fig. [8] we 
plot the clump radii versus the extinction-corrected clump ab- 
solute magnitudes. The typical clumps have radii of ~200 to 
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400 pc, while the DCOs have radii of 100-200 pc. Compar- 
ing the sizes obtained here with those obtained for the three 
DCOs using our more careful model simulations shown in 
Fig. |7] suggests that the more crude analysis may somewhat 
overestimate the true sizes of at least the DCOs. In any case, 
Fig. [8] shows an inverse correlation between clump luminos- 
ity and size, largely defined by the extreme properties of the 
DCOs. Their high luminosities (corresponding to very high 
stellar masses) and compact sizes imply very high densities. 
We will discuss the impUcations of this remarkable result later 
in the paper 

3.4. The relation between clumps and their host galaxies 

We have demonstrated above that massive star-forming 
clumps are an essential ingredient in our sample of local Ly- 
man Break Galaxy analogs, and have found that a subset con- 
tain a single extremely massive (up to several x K^Mq) com- 
pact (radius ^ IQ-^ pc) dominant central object (DCO). How 
do the properties of these compact structures relate to the 
overall properties of the host galaxies? 

The results are shown in the six panels of Fig. |9] In all pan- 
els we plot the extinction-corrected absolute S-band magni- 
tudes of the clumps measured within a fixed radius of 500 pc, 
in order to account for the effect that the angular sizes change 
by a factor of ^^2.5 over the redshift range z = 0.1 — 0.3. The 
results are as follows. 

Galaxy Concentration.- In panel (a) we plot clump AI b ver- 
sus the optical concentration parameter'"* of the host galaxy, 
C = 51og(^^). In general the more luminous clumps oc- 
cur in the LB As with larger concentrations. The six dominant 
central objects (filled squares) are the most extreme. Inspec- 
tion of Fig. [U shows that these objects are indeed systemati- 
cally found inside larger, optically faint disks compared to the 
sample on average. 

Galaxy Half-light radius- In panel (b)we plot M b versus the 
half-light radii for each galaxy measured in the UV images. 
This plot shows that the total UV light in galaxies hosting a 
DCO is dominated by the DCO itself, albeit with one notable 
exception (DCO object 080232 having Rjjv ~ 3 kpc). 

Stellar mass- In panel fcj we plot Mb versus the stellar mass 
of the hosts. The brighter clumps are predominantly found in 
the more massive objects in our sample. The DCOs are found 
in galaxies with masses A/, 2± 10'^'^^'^'^ Mq. 

Mid-IR luminosity- In panel (d)we show the relation between 
optical clump luminosity and the host galaxy luminosity at 
24/im. We find a correlation with L24^m indicating that the 
galaxies hosting the brighter clumps seen in the UV-optical 
are also the stronger IR sources. Even though selected on the 
basis of their high luminosity in the far-UV, the galaxies with 
DCOs have mid-IR luminosities similar to those of luminous 
infrared galaxies (LIRGs). 

To cal culate the concen tration we follow the procedures described in 
full detail in 'Lotz et al.' (2005') and Paper 1. First, we use SExtractor to make 
an object segmentation map and mask out neighboring objects. The image 
is background subtracted, and we calculate an initial Petrosian radius (rp 
with r; = 0.2) using the object center and (elliptical) shape information from 
SExtractor We then smooth the image by a = rp/5 and create a new seg- 
mentation map by selecting those pixels that have a surface brightness higher 
than the mean surface brightness at the Petrosian radius. We recalculate the 
object center by minimizing the second order moment of the flux, and then 
recalculate the Petrosian radius in the original image using this center. The 
total flux is defined as the flux within a radius of 1.5 X rp and C is then 
calculated in ckcular apertures containing 20 and 80% of the light. 



SFR.- These high mid-IR luminosities imply high star- 
formation rates (SFR^fQ+24), as shown in panel (e). The star 
formation rates are systematically larger in the LBAs with the 
more luminous clumps, and the six galaxies with DCOs have 
star formation rates ranging from 15 to nearly 100 A/q per 
year 

Comparing SFR indicators.- In panel (f) we plot Mb ver- 
sus the ratios of the two commonly used SFR indicators 
SFR/zq and SFR^fQ, ofcs+24- For the LBAs with the lower 
luminosity clumps, the majority show good agreement be- 
tween the two SFR indicators. The LBAs with the brighter 
clumps are sytematically displaced towards low values of 
SFRHa,a/SF^Ha.obs+24, with the galaxies having a DCO be- 
ing among the most extreme. To further investigate the dis- 
crepancy between SFR/ja o and SFRHa,ob+24 for the bright- 
est sources, we compare SFR^fa.o/SFR^fQ.obs+24 with the 
ratio SFR/fa o/SFRi^yv^.o- The result is shown in Fig. [TO] 
comparing the three SFR indicators. For the majority of 
sources all three indicators agree to within a factor of ~ 2 
to 3, but there are outliers for which the SFRs derived from 
the extinction-corrected Ha luminosity are significantly lower 
than either both or one of the other two SFRs. In other words, 
the ratio of the dust-corrected Ha luminosity relative to the 
luminosities of the mid-IR and (intrinsic) far-UV continua is 
systematically small in the galaxies with DCOs compared to 
normal star forming galaxies and the less extreme LBAs. We 
will refer to these offsets in Sections [4.11 and l4~2l below. 

4. DISCUSSION OF THE EMISSION LINE PROPERTIES OF LBAS 

The most significant results provided by our HST data are: 

1) Highly luminous, compact star-forming "clumps" 
appear to be the defining feature of the UV morphology of 
the LBAs. 

2) In the most extreme cases, the UV images are dominated 
by a single extremely bright and compact object surrounded 
by a disturbed envelope or disk. We have called these 
dominant central objects (DCOs). 

3) The star-formation rate derived from the extinction- 
corrected Ha emission-line luminosity tends to be system- 
atically smaller than that derived using the mid-IR and/or 
the extinction-corrected far-UV continuum luminosities (see 
panel/of Fig. |9]and Fig. [TOl i. More properly, the ratio of the 
Ha luminosity to either the mid-IR or far-UV continuum is 
smaller on-average than in typical star-forming galaxies. 

For our discussion below it is important to address first 
whether any of these properties could be explained by the 
contribution from a Type 1 (unobscured) or Type 2 (obscured) 
AGN. 

4.1. Are they AGN? 

4.1.1. Type 1 AGN 

Could the DCOs be Type 1 AGN in which we have a di- 
rect view of the accretion disk and broad emission-line re- 
gion? This would account for the exceptionally compact mor- 
phology as well as for the significant 24/im "excess". We 
note that the initial LBA selection process excluded any ob- 
vious Type I AGN in the first place based on the absence of 
obvious broad (several thousand km s^'^) Balmer emission- 
lines from the AGN Broad Line Region in the SDSS spectra 



021348 F150LP 



0.01 




[■WHM^ = (],U25" 
FWHM,„ = n.075" 



0.1 

Radius (arcsec) 




Overzier, R.A. et al. 
080844 F330W 




0.1 

Radius (arcsec) 



092159 F150LP 



10.00 



1.00 



0.10 



0.01 




Radius (arcsec) 



Fig . 7 . — Radial flux profiles of the three strongest, point-like sources identified 
ACS/SBC or HRC images compared to the measured profile of a model PSF (da. 
the size limits we can place on the nuclear sources, we simulated a range of 2D 
the model PSF. The results are shown using coloured lines. In all three cases, 
corresponding to 135-320 pc at 2 = 0.1 — 0.3). 



-17 r 



-18^ 

a -19 



-20: 

-21 \ 



□ □ 
□ □ 



I 



I 







100 200 300 400 



R 



50,clump 



(pc) 
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brightest clump in each galaxy estimated in Sect. 13.31 The six dominant 
central objects (DCOs) are indicated by filled squares. 



jHoopes et alj|2007l) . We can quantify this and exclude the 
possibility that some of our galaxies are of type Seyfert 1.8 or 
1 .9, with weak broad emission lines. In the event that possible 
broad emission was somehow being missed in the integrated 
SDSS spectrum, we have obtained integral field spectra using 
FLAMES on the VLT. In Fig. [TT]we show the preliminary 
results of the Ha+[NII] line profile as measured in the cen- 
tral ^1" region centered on the clumps in four LBAs host- 
ing some of the most luminous clumps (three of which are 
DCOs). The line profiles were fitted with a series of narrow 
and broad Gaussian components indicated by red and blue 
dashed lines, respectively. The fits show that each source has 
significant Une asymmetries on the blue side of the narrow 



in the HST images. Panels show the measured coun t rates (pluss es) in either the 
shed line) modeled using the Tiny Tim 3.0 software l'KrislF1995). To investigate 
Gaussian models having a FWHM of 0" 005-0'.' 125 and convolved them with 
the data is most consistent with an object no larger than 0'.'075 (FWHM; 



lines that are consistent with a relatively broad, blueshifted 
component (Overzier et al., in prep.). Such blue asymme- 
tries of width of a few hundred km s^^ in both the permit- 
ted and forbidden are not related to the much broader line 
profiles from Type 1 AGN, which are seen only for the per- 
mitted lines. Instead, the broad blueshifted wings seen on 
the emission-line profiles in our objects (Fig. fTll) are sign- 
posts of galactic winds and are commor ily observed in intense 
local and high redshift starbur sts (e.g. iH eckman et al."1990|; 
Lehnert & Heckman 1996; Lehnert et al...2009: Shapiro et al] 
2009h . 



4.1.2. Type 2 AGN 

Our spectra show clearly that the optical continuum must 
be dominated by the light of young stars (with no detectable 
contribution from an unobscured AGN). However, is it possi- 
ble that the DCOs are associated with the presence of an ob- 
scured (Type 2) AGN? In Fig. [I2]we show the log([NII]/Ha) 
versus the log([OIII]/H/?) diagnostic "BPT" diagram that is 
most commonly used for identifying the main source of ion- 
ization among samples of star forming galaxies and AGN 
([Baldwin et all [T98U iKauffmann et~aiTl200l iKewlev et alj 
I2006^ . Because it is based on ratios of bright emission lines 
that are observable out to z < 3, the diagram is widely 
used to characterize galaxies at both low and high redshift. 
It has been known for some time that intensely star-forming 
galaxies at high redshift such as Lyman breaks, "sBzK", 
and distant red galaxie s are often peculia r ly offset towards 
higher line ratios (e.g. 'Teplitz et al."2000"; 'Erb et al."2006a); 
l&ieketal. 2007; Lehnert etal. 2009; Hayashi et al. 200^ 
with respect to the mean SDSS star-forming population at low 
redshift (Fig. [121 dotted line). Similar offets, albeit perhaps 
slightly less extreme, have been observed in local "warm" 
IR-luminous galaxies (Kewlev et al.l I200"l1) and Wolf-Rayet 
galaxies (Brinchmann et al. 2008b ), and pEE P2 galaxies at 
z - 1.4 ( Shapley et alji2005bt iLiu et al.ll200 8). At present 
it is not known what physical mechanism s lie at the root of 
these o ffsets (see the detailed discussion in lBrinchm ann et aQ 
l2008ah . but it is important to note that our sample (filled 
circles) shows similar offsets in the BPT diagram (see also 
Hoopes et al. (2007) and Paper I). The BPT displacement is 
in the sense of enhancements in one or both of the two line 
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Fig. 9. — Correlations between the absolute B-band magnitudes (Mb) of the star-forming clumps and various properties of the host galaxies measured from 
the HST optical image, the SDSS spectrum and the Spitzer IR photometry. The six dominant central objects (DCOs) are indicated by filled squares. Panels show 
(a) Alg versus optical concentration index, (h) galaxy UV half-light radius, (c) stellar mass, (d) 24/^m luminosity, (e) total star formation rate, and (fl the ratio of 
the Ha SFR, SFR^f^ q. con'ected for dust using the Balmer lines to SFRj/jj o{,j,^^24fim ■ the SFR based on the attenuated Ha con'ected for dust using the 24/xm 
luminosity and Calzetti et al. ( 2007). We observe the following trends. The brightest clumps tend to have (a) larger host concentrations C, (b) smaller effective 
radii, (c) larger host stellar masses Af*, (d) larger host mid-IR luminosities, (e) larger host star formation rates, and (f' ) relatively low total SFRs based on Ha 
corrected for dust using the Balmer decrement compared to SFRs based on Ha corrected for dust based on the 24/^m emission from dust heated by young stars. 
The six dominant ce ntral objects (DCOs) are indicated by filled squares. The DCO having a relatively large UV half-light radius of ~3 kpc seen in panel b is 
object 080232. See ^^for details. 



ratios. If we focus on the locations of the host galaxies of 
the DCOs in Fig. [12] (filled squares), we see that they are 
more significantly displaced (in log[[NII]/Ha]) with respect 
to the locus of normal star-forming galaxies and the other 
LBAs. Typical galaxies in the SDSS with similar emission- 
line ratios are almost certainly composite objects in which the 
SDSS fiber encompasses both gas near the nucleus that is ex- 
cited by an AGN and n ormal regions of star formation in the 
surrounding galaxy (e.g.lKauffman n et al.l2003l : lKewley et al.l 
[2006; Kauffmann & Heckman 2009). Could the DCOs harbor 
a Type 2 AGN surrounded by a starburst? 

To answer this question we further investigate the BPT off- 
sets in Fig. [13] In panel faj we plot for each clump the di- 
mensionless parameter dp, defined as the relative offset its 
host galaxy has in the BPT diagram. This offset distance dp 
is measured perpendicular to the SDSS star-forming "ridge" 
(dotted line in Fig. \T2\ . This clearly demonstrates a relation 
between the offsets and the luminosity of the clumps. Panel 
(b) shows the offsets along log([NII]/HQ!) only, but we note 
that the correlation here is driven mainly by the well-known 
relation between mass and metallicity. In panel (c) we plot 
Mb versus the electron density (Ne) of the gas as derived 
from the density sensitive [SIIlAA6713A,67 3lA line doublet 
assuming a gas temperature of 10* K ( O sterbrock & Ferlandl 
120061) . The densities derived from the integrated SDSS spec- 
trum correlate with the Mb of the clumps over almost two or- 
ders of magnitude (the lowest and the highest being ~30 and 
^^2000 cm^'^) albeit with quite a large scatter The implied 
gas pressures are P/k - 10^ to 10^ K cm , several orders 
of magnitude higher than in normal star forming galaxies in 



the local universe. 

The Type 2 AGN interpretation would be consistent with 
the high electron densities seen in the most extreme cases 
(^ 10 cm^'^), since such densities are typical in the narrow 
emission-line region in AGN (Osterbrock & Ferland 2006). It 
would also be consistent with the blue-asymmetric emission- 
line profiles seen in Fig. [TT] whic h are characteristic of 
outflows in emission- Une AGN (e.g. iHeckman et akl [TosTl: 
[Whittle et aT][T988l) . The AGN hypothesis could also explain 
why the ratio of the 24 micron and extinction-corrected Ha 
luminosity is higher than in typical star forming galaxies (Fig. 
[Tol l, since AGN are characterized by strong emission by warm 
dust peaking in the mid-IR (e.g. Spinoglio et al. 2002). It 
would not however explain the result in Fig. [TO]that the ratio 
of the extinction-corrected far-UV and Ha luminosity is also 
higher on-average in the galaxies with DCOs than in typical 
star-forming galaxies, since in a Type 2 AGN the continuum 
from the ACjN itself will be very heavily obscured in the far- 
UV. 

To further test this starburst/ AGN composite hypothe- 
sis, we can use other emission-line diagnostic diagrams 
that are complementa r y to t he B PT diagram. Foll owing 
IVeilleux & O sterbrock] (il987h and iKewIev et ak (2006) we 
show in Figure [14] plots of [SII]/Ha vs. [OIII]/H/3 and 
[OI]/Ha vs. [OIII]/H/3. Bona fide starburst-AGN composite 
systems will lie in the region intermediate between the locus 
of pure star-forming galaxies and pure AGN in these plots. In- 
stead, we see that the galaxies with the DCOs tend to lie in a 
kind of "no man's land" displaced from the locus of star form- 
ing galaxies in a direction opposite to that of the AGN (they 
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were calculated using the recipes given in Sect. 12.41 Under reasonably 
standard assumptions for young starburst galaxies, all three SFR indicators 
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discuss the implications in Sections |4. 1 1 and The DCOs are marked by 
filled squares. 
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Fig. 11.— FLAMES/GIRAFFE/ARGUS spectra showing the Ha and 
[NlI]AA6548A,6584A line complex within a small aperture placed on the 
compact clumps in four of the sources. Sources 021348, 080844 and 210358 
host the most luminous and compact clumps seen in the LBA sample and all 
lie in the SF/AGN composite region of the BPT diagram making them can- 
didates for hosting a (weak) AGN. Source 005527 (bottom right panel) lies 
in the HII -region dominated part of the BPT diagram, eventhough it has one 
of the largest BPT offsets dp (see Fig. [12] and panel (a) in Fig. |9). The 
line profiles were fitted with a series of narrow and broad Gaussian compo- 
nents indicated by red and blue dashed lines, respectively. The fits show that 
each source has significant line asymmetries on the blue side of the lines that 
are consistent with a relatively broad, blueshifted component (Overzier et al. 
in preperation). Such blue asymmetries are typical for °alac tic outflows in 
starburst galaxies at both low and high redshift. See Sectr i4.ll for details. 

have weaker [SII] and [OI] emission than normal star form- 
ing galaxies). Their anomalous behavior is even more dra- 
matically illustrated in the right-most panel in Fig. [T4l where 
we plot [NII]/[SII] vs. [OIII]/H/3. As indicated in Figs. [12] 
and [14] (small grey circles), similar beh avior is seen in local 
IR-warm starbursts (' Kewlev et al.ll200l1) . 
Thus, the hypothesis that the compact clumps are found in 



starburst/ AGN composite systems does not explain the prop- 
erties of these objects. In future papers we will describe ob- 
servations with XMM-Newton, the Spitzer Infrared Spectro- 
graph, and the European VLBl network that are designed to 
search for more direct signatures of the presence of an AGN. 

4.2. Further Clues from the Emission Line Properties 

We have estimated the ages of the clumps and DCOs to 
be ^5-100 Myr assuming a model of an instantaneous burst 
of star formation (see Table |3]l. We have also shown that 
there exist striking correlations between the properties of the 
clumps (e.g., luminosity), the properties of their host galaxies 
(e.g., mass, radius, and concentration) and the physical prop- 
erties of the surrounding gas (e.g. ionization and density). 
The LBAs with DCOs are particularly extreme in these plots. 
Below we will discuss the most likely physical mechanisms 
that can explain these correlations. 

The results presented in Section 3 showed that there are 
significant correlations between the properties of the massive 
compact clumps and the physical properties of the surround- 
ing ionized gas as traced by the nebular emission-lines. In 
particular, the DCOs tend to be found in galaxies in which: 

1. The star-formation rate derived from the extinction- 
corrected Ha emission-line luminosity tends to be system- 
atically smaller than that derived using the mid-lR and/or 
the extinction-corrected far-UV continuum luminosities (see 
panel (f) in Fig. |9]and Fig. [TOl i. More properly, the ratio of 
the Ha luminosity to either the mid-lR or far-UV continuum 
is smaller on-average than in typical star-forming galaxies. 

2. The location of the galaxy in the diagnostic plot of 
[OIII]/H/3 vs. [NII]/Ha (BPT diagram. Fig. [HJ is more 
significantly displaced with respect to the locus of normal 
star-forming galaxies (see panels (a) and (b) in Fig. [TsT l. This 
displacement is in the sense of enhancements in one or both 
of the two above ratios. 

3. The electron density in the ionized gas is significantly 
higher than in the other galaxies (panel (c) in Fig. [I3] l. 

Occam's razor drives us to ask if the re is a single expla- 
nation for all of these results. In Section l^Tl we have already 
argued that the contribution from an obscured, low-luminosity 
ACjN could explain some but not all of the properties ob- 
served. Here we will consider two alternative possibilities: 
1) The DCOs are short-duration starbursts caught at a time 
after the peak in the O-star population, and 2) The DCOs in- 
clude "leaky" systems in which a large fraction of the Lyman 
continuum is escaping. 

4.2. 1 . The Aged Starburst Hypothesis 

The DCOs are extremely massive 10^ M^), compact 
(~ 10^ pc) objects. These high densities 10^ pc~^) 
imply very short dynamical times {R/v ^ 1 Myr). On simple 
grounds of causality, it is therefore plausible that these objects 
could have formed all their stars in a Myr, a timescale sig- 
nificantly less than the characteristic lifetimes of the O stars 
that produce nearly all the ionizing radiation in a starburst. In 
such a case, the ratio of the ionizing luminosity to the non- 
ionizing far-UV luminosity changes rapidly with time, and 
falls dramatically for ages greater than about 5-6 Myr (e.g. 
.Leitherer & Heckman .1995i) . The de-reddened UV/optical 
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Fig. 12. — BPT diagram of the galaxies in our sample (open/filled squares). Filled squai'es correspond to the LBAs hosting the DCOs. Compari son data 
shown are samples of star-forming galaxies at high redshift (green square: LBG cB58 iTepUtz et al. 2000); plusses: DEEP2 galaxies at 2 ~ 1 ( Shapley^eFal] 
l2005bl ;[Liu et al. 2008); red triangles: "BX" galaxies at 2 ~ 2 (Erb et al. 2006a); blue (red) upside-down triangles: star-forming (quiescent) X-selected galaxies 
at 2 ~ 2.5 (Kriek et al. 2007); dashed box region: star-forming "f,BzK" galaxies at 2 ~ 2 (Hayashi et al. 2009; Lehnert et al. 2009)) and low redshift (grey 
circles: nearby sample of warm IR-luminous galaxies from Kewley et al. (2001). Contours show the density distribution of emission line galaxies in the SDSS. 
Lines show the boundary relations between star-forming galaxies and AGN (Kauffmann et al. 2003, dashed line) and composite galaxies and AGN (Kewley et al] 
1200 6'. solid line), and the SDSS star- forming ridge line from lBrinchmann et al.l 12008a) . LBA 092159 (white cross) was found to contain a compact radio core 
using VLBl (in prep.). See Sect. |4. ll for details. 



colors of the massive compact clumps shown in Fig. |5](while 
uncertain) are consistent with those of such "aged" instanta- 
neous starbursts with ages of one-to-a-few tens-of-Myr (see 
Table|3]l. This is consistent with lBasu-Zvch et al.l (120071) who 
compared the SFRs derived from the UV, IR and radio. We 
should note however that at least in one case (210358) we find 
evidence for WR stars in the SDSS spectrum indicating that 
its SFH must include a much more recent (< 5 Myr) burst. 

If, however, the aged starburst hypothesis is correct for 
other objects, this would naturally explain the relative weak- 
ness of the Ha emission-line in these galaxies. Such a model 
could also explain the peculiar location of these objects in the 
BPT diagram (Fig. \T% . For an instantaneous burst in the 
age range between about ten and fifty Myr, the primary heat- 
ing/ionizing source for the associated gas would be the me- 
chanical energy supplied by su pernovae rathe r than O star ion- 
izing radiation (Leitherer & HeckmanPl995^ Models of such 
shock-heated gas (e.g. Dopita & Sutherland 1 995l;lAllen et al.l 
120081) show an enhancement in the [OIII]/H/3 and [NII]/Ha 
ratios, as we observe (Fig. [T2] |. 

The expected supernova rate would be ~ 4 (Mdump/^O^ 



Mq) per decade. Such a high supernovae rate occurring in 
a very small region also naturally accounts for the high elec- 
tron densities and thus the high gas pressures (P/k 10^ K 
cm^'') in these objects. Following Chev alier & Cle gg ( 1985), 
in the absence of significant radiative losses, the central pres- 
sure inside a starburst region of radius rioo (in units of 100 
pc) with a supernova rate of RsNe per decade is given by 
P/k = 1.4 X lO^RsNer^Qo K cm^^. iChevalier & Cleggl 
(1985) also demonstrate that this high-pressure gas will sub- 
sequently expand outward to form a galactic wind with a ram 
pressure that drops like r^^. This wind would be sufficient 
to produce the observed gas pressures over kpc-scales in the 
region surrounding the clump. The galactic wind will also 
accelerate emission-line clouds, leading to the broad blue- 
shifted wings seen on the emission-line profiles in our ob- 
jects (Fig. [TTI) . In fact, similar high gas pressures, broad 
blueshifted wings, and other signposts of galactic winds are 
common in intense local and high redshift starbursts (e.g. 
Lehnert & Heckmanll 1 996HLehnert et ai]|2009l: [Shapiro et alj 

This model fails however to explain the relative weakness 
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Fig. 13. — Correlations between the absolute B-band magnitudes (Mb) of the star-forming clumps and various emission line properties of the host galaxies 
measured from the SDSS spectrum. In each plot, the six dominant central objects (DCOs) are indicated by filled squares. Panels show (a) Mb versus BPT 
diagram offset dp, (h) log([NlI]/HQ:), and (c) elect ron d ensity A^e- The host galaxies of the brightest clumps tend to have larger BPT offsets dp, larger 
log([NlI]/Ha), and larger electron densities A^e. See ^4. 1 I for details. 
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Fig . 14.— In order to fu rther test the AGN hypothesis we show plots of [SII]/Ha vs. [OIII]/H/3 and [OI]/Hq vs. [OIII]/H/3 (following IVeilleux & Osterbrockl 
ll987tlKewlev et 311120061) . Bona fide starburst-AGN composite systems will he in the region intermediate between the locus of pure star-forming galaxies and 
pure AGN in these plots. Instead, we see some of the galaxies with the more massive DCOs lie in a kind of "no man's land" displaced from the locus of star 
forming galaxies in a direction opposite to that of the AGN (they have weaker [SII] and [OI] emission than normal star forming galaxies). Their anomalous 
behavior is even more dramatically illustrated in the rig ht-m ost panel where we plot [NII]/[SII] vs. [OIII]/H/3. Grey circles indicate the nearby warm IR-luminous 
sample from iKewlev et al. (200 L) . See the text in Sect. 14.1 I for details. 



of the [SII] and [Ol] lines seen in Fig. [14] Shock heating 
should lead to an enhan cement in the relative stre ngths of 
these lines (see Fig. 2 in Dopi ta & Sutherlandll 19951) . We are 
left with an otherwise attractive model that does not naturally 
explain all the properties of these objects. 

4.2.2. The Leaky Starburst Hypothesis 

An obvious way to lower the amount of observed Ha emis- 
sion would be a "leaky" interstellar medium in which a sig- 
nificant fraction of the ionizing radiation produced by O stars 
escapes the galaxy. This hypothesis has a plausible phys- 
ical basis: the feedback effect of concentrating the release 
of the kinetic energy from of-order 10^ supernovae (~ 10^* 
ergs) all within the exceptionally small volume of a DCO 
may create large holes in the surrounding interstellar medium 
through which ionizing radiation could then escape. Taken 
at face value, Fig. [TO] would then imply that the majority of 
the ionizing radiation is escaping in our two most extreme 
objects (0213 48 and 080844). It is interesting to note that 
IShapley et aH (|g006) find 2 out of 14 LBGs at z - 3 with a 
significant escape fraction (/esc — 50 — 100%) as measured 
from the ratio of the flux es below and above the Lyman limit 
in UV spectra. Similarly, l lwata et aTl (l2009h report the detec- 
tion of large /esc in 7 of 73 LBGs at z ~ 3. These fractions 



are similar to what is implied for our sample of LBAs under 
the leaky starburst hypothesis. 

The leaky starburst hypothesis can readily account for the 
high electron densities (gas pressures) and the broad blue 
wings on the emission-line profiles: these are still intense star- 
bursts and the same arguments as given above will apply. A 
leaky interstellar medium could also explain at least some as- 
pects of the peculiar emission-line ratios we observe. Models 
of photoionized gas under "matter-bounded" conditions (e.g. 
'Binette et al.' 1996: Casteflanos et al.ll2002t iGiammanco etal] 
2005; Brinchmann et al. 2008 show that the strengths of 
the [OI] and [SII] lines relative to Ha are decreased. This 
is because these lines are relatively strong in the warm par- 
tially neutral region lying beyond the edge of the Stromgren 
sphere, and this region is missing in matter-bounded condi- 
tions. These models show that the [OIII]/H/3 ratio will also be 
enhanced in matter-bounded clouds, as we observe (Fig. [T2] |. 
However, this model fails to explain the enhanced [NII]/Ha 
ratios seen for the DCOs in the BPT diagram. Models (B. 
Groves, private communication) show that matter-bounded 
clouds in starbursts should have [NII]/Ha ratios that are sim- 
ilar to or lower than ionization-bounded clouds. 

A more fundamental objection to the leaky starburst hy- 
pothesis is that is difficult to envisage a situation in which 
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a sizable fraction of the ionizing radiation escapes the galaxy, 
but the majority of the non-ionizing UV radiation is absorbed 
by dust inside the galaxy (as is required by the high mid-lR 
luminosities we see). For this reason, we regard the "aged 
starburst" hypothesis as the more likely one. 

In the future, we will use longslit and integral field spec- 
troscopy in order to spatially resolve and disentangle the line 
emission coming from the clumpy versus the diffuse regions 
seen in the HST images. 

5. DISCUSSION OF THE STRUCTURAL PROPERTIES OF LBAS 

The peculiar emission line properties discussed above show 
the presence of significant feedback processes between the 
massive star-forming regions in LBAs and their surrounding 
ISM. Here, we will focus on the remarkable stellar properties 
of the massive clumps and the DCOs in particular, by asking 
the question to what kind of structural components of galaxies 
do they correspond? 

The DCOs are only marginally resolved in the HST UV im- 
ages, with implied radii of no greater than '^100 pc (see Fig. 
I?]). The fainter clumps are typically several times larger than 
this. The masses of the DCOs estimated in Section 3 range 
from a few x 10^ to a few x 10° M©, while the more typical 
clumps have masses in the range of a few x 10^ to a few x 10* 
Mq. This implies very high stellar mass surface densities. For 
the typical clumps values for the effective surface mass den- 
sity (Ee ~ M*,ciurnp/i2nRl)) range from - lO* to 10° M© 
kpc~2^ while for the DCOs Eg - lO^-^ to 10" Mq kpc'^. 
What then will the clumps and DCOs evolve into? Can we 
identify such evolved counterparts in local galaxies? 

5.1. Super Star Clusters and Central Massive Objects? 

Massive clusters are observed in and around the nu- 
clei of many nearby late- and early-ty pe galaxies, mergers 
and merger remnants, and AGN (e .g. Carollo et al. 1997; 
Geha et all 120021: iBekeretal] 120041: iMaraston et al. 2004; 
Walcher et al.ll2005l:ISeth et alJ2008h . The largest stellar clus- 
ters typically have sizes of only a few pc (the largest being 
« 10 pc) and masses rarely exceeding ~ 10^ Mq. Their struc- 
tural properties are thus similar to those of massive globular 
clusters in the Milky Way and local galaxies, and also to the 
young "super star cl usters" (SSCs) found in starburst galaxies 
(IMeurer et al]IT995h . 

The most luminous star clusters in normal and starburst 
galaxies he along a relation between the absolute optical mag- 
nitude^ oi^ the cluster and the global SFR of the host galaxy 
(|Larsen 2002), a relation that is mainly driven by the rela- 
tion between the luminosity of the brightest clus ter and the 
total number of clusters formed in a galaxy ( Wh itmore et alJ 
I2007h . The correlation implies that for galaxies having SFRs 
as high as LBAs even the brightest clusters should only have 
My ~ -16 to -14 (My = -1.871og(SFR) - 12.14, see 
IWeidner et all |2004l) compared to -21.5 < My < -18 
found in Fig. |5] This discrepancy could be partly reconciled 
given that the LBA clusters are systematically younger than 
the brightest clusters observed in normal galaxies, since a star 
cluster fades by ^ 3 mag in the first few 100 Myr. However, 
this scenario still largely falls short for any of the luminous 
clumps/DCOs having -21.5<Afy<-19. 

The most massive cluster-like object known to date, object 
W3 in the merger remnant NGC 7252, has a mass of (8 ± 
2)10^ M© and a radius of 17.5±1.8pc ( Maraston et alj2004l) . 
Interestingly, its high mass and small size are very much at 



odds with any kind of stellar cluster found to date. Curiously, 
the properties of object W3 are perh aps more similar to that 
of ultra-compact dwarf galaxies (e.g.' Mieske et alj|2004l) . If 
we extrapolate from its current luminosity at an age of ~400 
Myr to t hat expected at 10 Myr, we find M^°*^^'' ~ -19.2 
( lBastianll2008h . still over 2 magnitudes fainter than some of 
the DCOs. 

Although the stellar mass surface densities that we derive 
(~ 10*-^ - 10" Mq kpc^^) are comparable to that of typical 
massive clusters mentioned above, the sizes are about an order 
of magnitude larger (^100 pc) and the masses are up to two 
orders of magnitude larger (< 10° Mq). Thus, the massive 
clumps and DCOs identified in the LBA sample appear as sig- 
nificantly scaled-up versions of known types of massive star 
clusters, and are not simply an extension of the cluster lumi- 
nosity function as derived from nearby star-forming galaxies. 

Furthermore, nuclear super star clusters and supermas- 
sive black holes are jointly referred to as central mas- 
sive objects (CMOs). CMOs follow a tight correla- 
tion with the total mass of the galaxy over many orders 
of magnitude in stellar mass. The central star clusters 
in low luminosity early-types and spiral bulges, as well 
as the masses of the central supermassive blackholes in 
L > L^, early-types all lie on t he relation Mcmo ~ 
0.002A/oa; with small scatter (e.g. 'Magorrian et al.' '19981; 
Marconi" & Hunt 2003; Cote et al. 2006; Ferrarese et al. 200|; 
IWehner & Harrisi |200 6; Rossa et alj|200 6). Simulations of 
both isolated and me rging galaxies reproduce the observed 
correlation well (e.g. iLi et al.l 120071) . As we discuss below 
in more detail, the clumps and DCOs typically represent sev- 
eral percent of the total stellar mass of the surrounding galaxy. 
These mass fractions are therefore over an order of magnitude 
higher than in the CMOs. 

Could the massive clumps or DCOs consist of multiple 
SSCs? If so, how long would it take this ensemble to merge 
into a single object? To address this, we assume that a 
clump of mass 10* Mq actually consists of 10-100 SSCs 
each having a mass of lO^"'^ Mq distributed over a re- 
gion with a 100 pc radius within which the velocity disper- 
sion (Jciump « \/GMciump/5Re « 30 km s"^ Then the 
timescale over which a single SSC migrates inward as a result 
of dynamical friction is on the order of only ^10-20 Myr 
(|Binnev & Tremaine 2008), comparable to the typical ages 
we have derived for the clumps. This suggests that a given 
massive clump (or DCO) is likely to have evolved into a sin- 
gle very massive object, even if it started out as an ensemble 
of SSCs. 

We conclude that the clumps and DCOs are very different 
structurally from typical nuclear (super) star clusters. 

5.2. Extra light and central "cusps" in elliptical progenitors? 

For the discussion below it is important to point out that 
the DCOs (and, more generally, the most luminous clumps) 
tend to appear in the more massive galaxies in our LBA sam- 
ple (see panel c in Fig. |9]l. The DCOs also each sit in the 
middle of extended, perturbed disk-like structures (hence the 
correlation between clump luminosity and the optical concen- 
tration index, panel a of Fig. |9j. Another characteristic (by 
definition) is that the DCOs are single objects rather than the 
multiple (fainter) clumps as seen in many of the other objects. 

In Fig. [15] we plot the clump mass (Af,) and effective ra- 
dius (_Re) versus the effective mass density (Eg) for the bright- 
est clump in each LBA, and compare with a collection of 
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Fig. 15. — Mass (M*) and effective radius (Re) versus the effective mass surface density (Se) of the brightest clump identified in each galaxy (squares). The 
compaiison data shown was extracted from a compilation of Hterature results presented in Fig. 45 in H opkins et ah 1,2009). showing globula r clusters (crosses), 
nuclear star clusters (triangles), and the extra-light components in early-type galaxies from the samples of lKormendv et alj j2008, stars) and lLauer et alj (20071 
circles) as measured by Hopkins et al. 1 2009). Objects that we identify as DCOs based on their location in the M«, Re, Se plane in combination with the large 
total stellar masses and concentrations of their host galaxies are indicated with filled squares. We suggest that DCOs could be the progenitor starbursts of the 
central "extra light" component in local cuspy elliptical galaxies. See Sect. l5.2l for details. 



results from the literature. The comparison data shown was 
extracted from a compilation of literature results presented 
in Fig. 45 in Hopkins et al. (2009), showing globular clus- 
ters (crosses), nuclear star clusters (triangles), and the extra- 
light components in ea rly-type g alaxies from the sa mples of 
iKormendy et alj ( 120081 stars) and lLauer et"!!] (120071 circles). 
The high masses of ^ 10^ Mq derived for the DCOs com- 
bined with the strong constraints on their small effective radii 
(^ 100 pc as evidenced by the HST diffraction spikes in Fig. 
m yield the highest mass surface densities found in our sam- 
pl e, Sm fti 10^° - 10" Mq kpc-2 (see TableH. 

iHopkins et al.l (120091) show that this region of parameter 
space in the M* and Re versus Sg planes is almost exclu- 
sively occupied by the so-called central excess or "extra light" 
component seen in local el liptical g alaxies (star s and circles in 
Figj31 ICrane et al.lll993HLauer etld]fT995l:lKormendy et al.1 
|2008[ and references therein). This extra light is characterized 
as an inner (<1 kpc) component that rises more steeply than 
expected from the extrapolated Sersic power-law that fits the 
outer profile of ellipticals well. The so-called "cusps" have 
effective radii of 50 — 500 pc and show a trend of decreas- 
ing density with increasing mass and in that sense they lie on 
the same sequence formed by elliptical galaxies. The cuspy 
cores are furthermore typical for relatively low-mass (< L,) 
ellipticals, while the y are relatively u ncommon for the most 
massive ellipticals. Kormendv et al.l (|^08) suggest that the 
core-cusp dichotomy in local ellipticals may correspond to 
their formation history: the extra light or cusp components 
are very likely the remnant signature of a massive central star- 
burst in a dissipative ("wet") merger, while the extended cores 
in core-dominated ellipticals are consistent with being the re- 
laxed descendants of dissipationless ("dry") mergers. 

While none of the clumps bear any structural similarities 
with both globular and nuclear star clusters (crosses and tri- 
angles), the DCOs in the M*, Re, Eg plane appear structurally 
similar to the extra-lig ht light components in early-type galax- 
ies (stars and circles). IHopkins et al.l (120091) investigate in de- 



tail how the central mass associated with the extra light com- 
ponent is related to the remnant of a massive central starburst 
fueled by dissipative merging, and use extensive simulations 
to show that such a mechanism successfully reproduces the 
observed profiles of elliptical galaxies ( Hopkins et al.ll2008l 
120091) . In particular, they find a strong correlation between 
the galaxy stellar mass and the fraction of the total stellar mass 
that is in the extra light component (f extra), as well as a one- 
to-one correlation between the extra light mass fraction and 
the fraction of the total stellar mass that formed in the dissipa- 
tive starburst if starburst), with roughly a factor of ~ 2 intrin- 
sic scatter in f starburst and a factor of ~ 4 scatter in /extra 
at each M*. For stellar masses in the range of 10^° to 10^^ 
Mq the fstarburst{f extra) ranges from a few to 50%, consis- 
te nt with the extra-light fractions measured in real ellipticals 
bv lHopkins eTal] ([2009h . 

We can compare this with the ratio of the central clump 
mass to the total stellar mass we found for the DCOs. Due to 
the uncertainties in the mass estimates involved, we estimate 
the DCO mass by taking the average of the values derived 
using single burst and continuous SFHs. The result is shown 
in Fig. [16] We find mass fractions typically of a few percent 
(J DCO ^ 2 ~ 4%), object 021348 being a notable exception 
Udco ^ 23%, see also Table |4|. Except for object 021348, 
which agrees well with the cuspy cores, the mass fraction of 
the DCOs appears to be relatively low compared to the early- 
type cusps plotted in Fig. [T6](stars and circles). However, the 
mass fractions we find are still significant and not inconsistent 
with the data. 

Alternatively, we should note that the special region of pa- 
rameter space in M*, Re and Sg occupied by both the DCOs 
and the cuspy cores appears to represent a particular tran- 
sition region between massive, cluster-sized objects on one 
hand and entire galaxies on the other. While low-mass star 
clusters follow the scaling relations of globular clusters for 
masses < 10^ A/q, a wide variety of higher mass objects 
ranging from the most massive nuclear star clusters, dwarf- 
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Fig. 16. — Fraction of the stellar mass in the DCOs (see Fig. |4) versus 
the total stellar mass of the six objects. The DCOs contain 2-23% of the 
total stellar mass, suggesting that they could be the hypothesized progenitor 
starbursts related to the central "extra light" component observed in local 
cuspy elliptical galaxies. This is also consistent with their total stellar masses 
of 10^" < AI^/Mq < 10^^ . T he comparison data shown (stars and circles, 
taken from Fig. 20 in lHopkins et al. 1 2009)) coiTespond to the same extra- 
light components as plotted in Fig. [15] See Sect. l5.2l for details. 

galaxy transition objects, ultra-compact dwarf galaxies, dwarf 
elliptical nuclei, and the early-type cusps themselves all ap- 
pear to l ie on the same scaling relations for spheroidal galax- 
ies (e.g. iKlssler-Patig et al]|2006t iHopkins et al.ll2009l) . Al- 
though this overlap may be partially coincidental, perhaps it 
can be explained if these objects are all the product of some 
generic dissipational process. 

In contrast with the DCOs, the much less extreme clumps 
indicated in Fig. [15] are typically found in LBAs contain- 
ing multiple clumps distributed within their central regions 
(see Figs. [T] and |3]l. For a clump of mass Mdump ^ 10* 
Mq moving within a disk of mass M^isk ~ 10^" Mq and 
a disk radius of Rdisk ^1-5 kpc (the average optical half- 
hght radius) we find adisk ~ y'GMdisk/iARdisk w 80 km 
s~^ and a dynamical friction time of ^^100-200 Myr. This 
is several times longer than the ages that we have derived 
for the individual clumps. Thus it is plausible that at this 
relatively early dynamical stage we usually see several such 
clumps within the effective radius of the galaxy. The diffuse 
light seen in the HST images may build up as the clumps lose 
stars due to tidal stripping while they interact with each other 
and with the disk and slowly migrate inward. This process 
may eventually lead to the formation of an inner bulge as well 
as the Sersic profile typical of disk galaxies (Noguchi 1999; 
Immeli et al."2004'; 'Carollo et al."1997t iBournaud et aUlOOli 
2008; Elmegreen et al. 2008a b, 200^ 

5.3. Relation to Mergers and IR- Luminous Galaxies 

Although most LBAs were unresolved in the original SDSS 
images, some showed evidence for close companions or faint 
extended emission. One of the major results from the HST 



images is that most of them show very complex morpholo- 
gies at scales well below the SDSS seeing or sensitivity. We 
interpret the large array of morphologies in Fig. [T] as ev- 
idence that most LBAs are undergoing or have undergone 
an interaction or merger'^. We also see a small number of 
objects that do not show unambiguous evidence of merging 
such as tidal tails or companions, but that are nonetheless of- 
ten highly peculiar as evidenced by ring or chain-like mor- 
phologies (e.g. object 102613 and objects 082001, 082550, 
083803). These observations suggest that the starbursts in 
LBAs and the formation of the massive stellar clumps are 
somehow linked to t he mergers (e. g. Larson & Tinslev 191% 
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Owers et al.ll2007t iLi et al.ll2008l) . The LBAs must also be 
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sufficiently gas-rich in order for any starburst trigger to be ef- 
fective. Simulations show that the inflow of gas ensures that 
the pressure in the interstellar medium becomes larger than 
the internal pressure of giant molecular clouds, and that the 
clouds can collapse and form massi ve clusters be fore they are 
disrupted by supernova explosions (: Bekkill2()04l) . 

It is important to note that the spectral energy distributions 
of LBAs are very different from typical starburst galaxies in 
the local Universe. These are predominantly associated with 
the highly dust-obscured LlRGs (10" < Lir < 10^^ Lq) 
and ULlRGs {Lir > lO^^ Lq). At least ULlRGs are known 
to be strongly associated with the different phases of galaxy- 
galaxy mergers, with about equal fractions showing dual nu- 
clei (pre-merger) and showing indication of post-merging (af- 
ter the separate nuclei have coalesced and a period of violent 
relaxation starts forming an ex tended spheroid a l component) 
(e.g. [G enzel et al. 2001; Tac coni et al.l 120021: lOasyra et al.l 
l2006allbi : lArribas et al.ii200 8). These mergers are found to be 
limited to galaxy pairs having near-equal mass ratios ( 1 .5 : 1 on 
average, with a maximum of 3:1). The less luminous LlRGs 
are typically associated with much less violent events that are 
less effective at triggering a major starburst phase as observed 
in ULlRGs, perhaps as a result of smaller gas fractions, or- 
bital parameters, and l arger pair mass ratio s (or a combination 
thereof). As shown in Arribas et al.l (12008 '), about one third of 
LlRGs appear in single isolated (spiral) galaxies. While none 
of the LBAs qualify as ULlRGs, there is some overlap with 
the LlRGs since -30-50% of LBAs have Lir > 10" L©. 
As such, LBAs likely represent a population of relatively un- 
obscured mergers or merger remnants that share some char- 
acteristics with LlRGs. The number density of LBAs is much 
lower than that of LlRGs. If the two populations are somehow 
connected, the LBAs could represent a very short evolution- 
ary phase during which the starburst can be seen relatively un- 
obscured. Alternatively, or in addition, LBAs might be seen 
at a special viewing angle, perhaps along the axis where the 
starburst has blown a hole in the ISM and is leaking Lyman 
continuum photons. 

5.4. Implications for Black Hole Formation 

We have shown that the LBAs are consistent with galaxies 
that are in the very early stages of forming their central bulge 
and/or core excess. Local galaxies have tight correlations be- 
tween the mass of their central black holes and that of their 
bulge and total mass. It is thus interesting to ask whether there 
is any evidence that a large black hole mass is accumulating 

" Here we use the term "merger" in the widest sense of the word, i.e. 
it is not limited to "major mergers" but includes minor mergers and other 
mechanisms that may deliver a supply of gas-rich material. 
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in the LBAs as well. The DCOs are particularly intriguing in 
this regard: such very massive and compact objects appear to 
be the ideal breeding gounds for a supermassive black hole. 
As we have no means of directly detecting the gravitational 
signature of a black hole inside the massive clumps, the only 
way to infer the presence of a massive black hole would be 
through the spectral signature of an actively accreting AGN. 
As discussed in Sect. 14 . 1 1 we did not find any direct evidence 
for energetically significant AGNs in the DCOs. 

It has long been suggested that the mass loss from stars and 
supernovae in a very massive compact cluster at the center 
of a galaxy may play a major rol e in the f ueling of a su- 
permassive black holes (Norman & Scovillei ,1988 ). While 
the stellar winds from O and WR stars and the ejecta from 
core-collapse supernovae may in principle provide a signif- 
icant amount of material available for accretion during the 
early stage of a massive central starburst, these are not ex- 
pected to be viable fueling mechanisms because the speed 
of the outflowing stellar ejecta i s far in excess of the es - 
cape velocity fro m the star cluster ( Norman & Scovillell988l) . 
iNorman & Scov ille ( 1988) argue that the high mass-loss rates 
and the low mass-loss speeds of the winds from evolved AGB 
stars may offer a steady flow of material that will be retained 
within the potential well of the star cluster. It can then cool 
and (some fraction) can be accreted by the supermassive black 
hole. Thus, the growth phase of the supermassive black hole 
(a powerful AGN) is expected to occur only after the time 
that the mechanical feeedback energy input from core col- 
lapse supernovae has subsided and the progenitors of AGB 
stars have reached the end of their main sequence lifetimes. 
This will occur at ^50 Myr after the onset of the starburst 
dLeitherer & Heckmanll 1995b . and this age constraint is con- 
sistent with the observed ages of the nuclear starbur sts ob- 
served in a sample of local AGN (iDavies et al.ll200'7h . Un- 
der the assumption of an instantaneous burst, the ages of the 
DCOs shown in Fig. |4] and listed in Table |4] are generally 
much lower (5-29 Myr) than 50 Myr. Thus, it is perhaps not 
unexpected that under such a scenario the DCOs are currently 
not capable of efficiently fuelling a supermassive black hole. 

We also note that in th e "migrating-clump model" advo- 
cated by Bourn aud et al.l (l2007l 120081) and Elmegree n et al.l 
dSoOSb), supermassive nuclear black holes may form under 
the assumption that each of the clumps gave rise to an in- 
termediate mass black hole prior to its migration toward the 
center where they can coalesce. Such a mechanism may be 
able to explain the presence of seed black holes of substantial 
mass in the centers of gas-rich galaxies dominated by mas- 
sive, turbulence-induced star-forming clumps. Both mecha- 
nisms described above may therefore predict significant delay 
times between the onset of a starburst and that of AGN. 



6. SUMMARY AND CONCLUDING REMARKS 

We have presented the results of the analysis of a sample 
of 30 low-redshift (z < 0.3) galaxies selected from the union 
of the GALEX and SDSS surveys to have the high fai--UV lu- 
minosities and small sizes characteristic of the high-r edshift 
Lyma n Break Galaxies jHeckman et al. 2005; Hoop es et al.l 
120071: iBasu-Zvch eTaTI 12007; Overzier et al.. ,2008.) " These 
"Lyman Break Analogs" (LBAs) offer the opportunity to in- 
vestigate the physical processes occurring in Lyman Break 
Galaxies (LBGs) in considerably more detail than is possi- 



ble at high-redshift . In this paper we reported on the inter- 
nal structures of the LBA's using HST optical and UV images. 
We complemented these data with mid/far-lR photometry and 
optical spectroscopy from the SDSS and the VLT. 

We have shown that the morphologies of the LBAs are 
dominated by massive compact clumps of young stars. Our 
most unexpected and potentially important result is that we 
have found six galaxies in which the galaxy UV morphology 
is dominated by a single very compact object which is unre- 
solved or only marginally resolved in the HST images. These 
dominant compact objects ("DCOs") have estimated masses 
of few xlO^ to few xlG^ M© and radii of-order 100 pc or 
less. In the other LBAs we instead find multiple clumps per 
galaxy, the brightest of which have estimated masses of typi- 
cally a few X 10* il/0 and radii of several hundred pc. 

We show that the galaxies containing the DCOs differ from 
the other LBAs in many respects: 

1. Their host galaxies tend to be more massive than the 
other LBAs (median masses of 3 x lO^*' vs. ^ 6 x 10^ 
Mq respectively). 

2. They have systematically higher central concentrations 
of light than the other LBAs (the DCO is typically located 
near the center of a larger disk or envelope). 

3. The DCOs are found in galaxies with high mid-infrared 
luminosities: L24^m ~ 10^" '^ to ~ 10^^'^ Lq, roughly an 
order-of-magnitude more luminous than the other LBAs and 
comparable to luminous infrared galaxies. 

4. The galaxies with a DCO are more strongly displaced 
from the locus of normal star forming galaxies in the standard 
BPT diagnostic emission-line ratio diagram. This displace- 
ment is in the sense of having larger [OIIIJ/H/3 and/or larger 
[NII]/Ha ratios. 

5. The DCOs tend to occur in galaxies with higher gas den- 
sities and pressures as traced by the ionized gas {P/k ~ 10^ 
cm-3 K). 

6. The DCOs are found in galaxies in which the 
star-formation rate inferred from the dust-corrected Ha 
luminosity tends to be smaller than that derived from the 
mid-lR luminosity and/or the extinction-corrected far-UV 
luminosity. In the most extreme cases, this discrepancy 
reaches roughly an order-of-magnitude. 

We have shown that the DCOs are not Type 1 (unob- 
scured) AGN. We find no evidence for the presence of a 
very broad (several thousand km/s) Ha emission-line, with 
the most stringent constraints provided by our VLT GI- 
RAFFE/ARGUS IFU spectra. Instead, these data show blue- 
asymmetric wings on the Ha and [Nil] line profiles indicative 

'^IHeckman et al.' '2005) coined the term "living fossils" in relation to the 
Lyman break analogs sample then discovered. Indeed, an apt analogy in the 
field of evolutionary biology is presented by the Coelacanth. A modem-day 
analog of this species of fish, known from the fossil record and believed to be 
extinct since the end of the Cretaceous period, was discovered in 1938. The 
find illustrates the use of analogs in the study of cosmology. Even though the 
local analogs will differ from their distant look-a-likes to some level of detail, 
and even though their local environment has probably changed drastically, 
they offer an entirely new way of investigating the fossil record alongside 
with other, more traditional fines of research. 
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of a dusty outflow of ionized gas. While we can not entirely 
exclude the possibility that an energetically significant Type 
2 (obscured) AGN is present in some or all of the DCOs, we 
regard this possibility as unlikely. In particular, the DCOs 
are associated with galaxies having unusually weak [SII] and 
[OI] emission-lines compared to normal star forming galax- 
ies. This is the opposite of what is seen in starburst-AGN 
composite systems (e.g. ,Kewlev et alj 120061) . but similar to 
what is seen in low redshift IR-warm starbursts jKewlev et alj 
l200lh . 

The continuum of the DCOs is dominated by the light of 
young stars: these are clearly strong highly compact star- 
bursts. This is consistent with both the high gas pressures 
and emission-line profiles: both are commonly observed in 
starbursts and are the signposts of supernova-driven galactic 
winds (e.g. Lehnert & Heckman 1996). We have considered 
two specific scenarios under which the starbursts associated 
with the DCOs might be unusual. They may be evolved short- 
duration starbursts in which the rate of photoionization is rel- 
atively low because few O stars are present. They may also be 
"leaky" starbursts in which a significant fraction of the ioniz- 
ing radiation escapes (matter bounded conditions). The latter 
scenario seems less likely, since it is difficult to envision cir- 
cumstances in which most of the ionizing radiation escapes 
the galaxy, but most of the non-ionizing UV radiation is ab- 
sorbed by dust and reprocessed in the mid/far IR. 

We have compared the properties of the DCOs and other 
massive clumps to those of the nuclear star clusters typically 
seen in different kinds of nearby galaxies. The DCOs are 
highly massive and dense objects that do not appear to be re- 
lated to single, young star clusters which are significantly less 
massive (< 10^ M0) and smaller (^10 pc) than the struc- 
tures we see in the LB As. However, the DCO structures are 
quite similar to predictions for the merger-driven starbursts 
that are believed to be the progenitors of the extra light or cen- 
tral cusp s seen generically in lower- mass (< M^,) el liptical 
galaxies jKormendv et al.ll2008t iHopkins et al.1 120091) . This 
would be consistent with our optical images of the LBAs, 
most of which strongly suggest on-going mergers. However, 
at present it is unclear whether the DCO hosts could passively 
evolve to form a fuU-grown elUptical without additional merg- 
ing. 

More generally, the properties of the LBAs are quite con- 
sistent with the idea that large Jeans instabilities in a gas- 
rich disk will lead to the formation of highly massive star 
forming clumps. At high redshift, it is not yet clear what 
the main source of this gas is: it could be replenished con- 
tinuously, or come in the form of discrete accretion events 



of small galaxies or clouds (or, most likely, a combina- 
tion thereof). These processes have been proposed to ex- 
plain the properties of high-z star-forming galaxies (e.g., see 
Overzier et al. 2008; Basu-Zvch et al. 2009; Elmegreen et al] 
l2009t IForster Schreiber et all |2009|; iLaw et al.i i2007l 120091 
and references therein). In any event, it is expected that the 
massive stellar clumps will aggregate at the center of a galaxy 
via dynamical frictio n to ultimately form a b ulge and super- 
massive black hole (lElmegreen et alJ [20091 and references 
therein). The LBAs offer us the opportunity to study such 
processes in more detail, and compare samples across differ- 
ent redshifts. 

We speculated about why we do not see clear evidence for 
the formation or rapid growth of a supermassive black hole 
in the DCOs (which appear to be the ideal black hole nurs- 
eries). One intriguing possibility is that these objects are sim- 
ply too_^oung_Jor_a_black hole to be forming yet. Follow- 
ing |NormmAlcoyille (1988) and Davies et al. (2007), the 
growth of a supermassive black hole inside a very massive 
compact star cluster may not commence until the cluster is old 
enough for the core-collapse supernovae phase to end (cluster 
age > 50 Myr). Subsequent to this, gas shed by evolved stars 
will be retained in the cluster's potential well rather than being 
blown out. This gas can then cool, flow inward, and be used to 
grow the black hole. If correct, we should be able to identify a 
population of the descendants of the young DCOs in a phase 
as powerful AGN. This scenario might also account for the 
rarity of strong AG N in the LBG population (e.g. Steide Tetan 
12002; Ouchi et al. 2008). 

Follow-up UV spectroscopy using the HST Cosmic Origins 
Spectrograph has been allocated for selected targets in Cycle 
17. The UV spectra will allow us to study many of the fea- 
tures discussed in this paper that are relevant to LBAs and to 
starbursts in general. These include the role of galactic out- 
flows and low-luminosity AGN, the escape of ionizing radia- 
tion, the strength of Lya emission, stellar abundances and the 
initial mass function. 
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Properties of the Lyman break analogs sample. 
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Optical half-light radius. 
Optical concentration index. 

Gas velocity derived from the Ha omission line width. 
Dynamical mass M^^.j^ = 3.4Kecr^/G. 
Dynamical time t^y^ = R^/cr^. 

Oxygen abundance in units of 12 + log{0 / H) estimated usina the "OSNT" estimator fron JPettini & Pagel[2Q04 l. 



Local Analogs to Lyman Break Galaxies 



21 



TABLE 2 

Emission line ratios and electron density. 
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TABLE 3 
Properties of the clumps. 
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t ID number of the host galaxy of each clump. 

t Magnitudes and colors correspond to those measured in the filters F850LP and F330W-F850LP, 
instead of, F606W and F150LP-F606W, respectively. 

'•^ Absolute (optical) clump magnitude from Fig.|5] M = nicor — 51og ^^i^) + 2.51og(l + z), 
where is the luminosity distance, rricor is the apparent magnitude corrected for dust using 
mcor ^m-E{B-V),k'{X) and assuming £(B - l^). = 0A4:E{B~V)gas (Calzetti 20011). 
Dust-corrected clump color from Fig. [5] 

Clump age and stellar mass under the assumption of the instantaneous star formation history 
shown in Fig. |5](red tracks). 

Clump age, stellar mass and star formation rate under the assumption of the continuous star 
formation history shown in Fig.[5](blue tracks). 
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TABLE 4 

Properties of the massive, unresolved clumps ("DCOs") discussed in Sect.[51 
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11 


10.3 


0.23 
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0.46 


0.29 


9.11 


202 


9.70 


74 
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50 


10.6 


0.03 


080232 


0.20 


0.17 


9.18 


156 


10.00 


91 


2 


153 


10.5 


0.04 


080844 


0.77 


0.30 


8.21 


95 


9.46 


38 


2 


56 


10.0 


0.02 


005439 


0.43 


0.23 


8.77 


190 


9.41 


52 


4 


70 


10.5 


0.02 



t ID number of the host galaxy of the DCO. 

Fractions of DCO-to-total light meaured in the UV (fuv) and optical if opt)- A circular aperture of 0'.'2 radius was used for the DCO, 
and an aperture of 10" radius was used to estimate the total light. Companion galaxies were masked out and a large object-free annulus near 
the source was used for background subtraction. For object 080844 the UV/optical light fractions correspond to those measured in the filters 
F330W/F850LP, instead of F150LP/F606W. The total optical light for 080844 was measured within an aperture of 5". 
Average of the single burst and continuous SFHs clump masses of Table[3] 
Physical optical half-light radius measured within the 0'.'4 diameter aperture used in Sect. 13.31 
= Stellar mass surface density Ea/ ~ {Af,,c!timp)/27r_R^ 



clump ■ 



G(A/.. clump) 



' Clump velocity dispersion assuming spherical geometry (Jv,ciump = y 
^ Dynamical time tdyn = Rciu,np/iyv,ci ump 

H Dynamical friction time = (p^) {^^)\^^^) , where A = 

' Galaxy stellar mass estimated from the average of the SED and dynamical masses, 
■i Stellar mass fraction of the clumps. 



